Electronic structure of unidirectional superlattices in crossed electric and magnetic 

fields and related terahertz oscillations 
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We have studied Bloch electrons in a perfect unidirectional superlattice subject to crossed electric 
and magnetic fields, where the magnetic field is oriented "in-plane", i.e. in parallel to the sample 
plane. Two orientation of the electric field are considered. It is shown that the magnetic field 
suppresses the intersubband tunneling of the Zener type, but does not change the frequency of 
Bloch oscillations, if the electric field is oriented perpendicularly to both the sample plane and 
the magnetic field. The electric field applied in-plane (but perpendicularly to the magnetic field) 
yields the step-like electron energy spectrum, corresponding to the magnetic-field-tunable oscillations 
alternative to the Bloch ones. 

PACS numbers: 78.45+h, 73.21. Cd, 73.40.-c,78.67.De 
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I. INTRODUCTION 

Semiconductor structures are considered as a perspec- 
tive source of persistent terahertz radiatiouji In semicon- 
ductor superlattices the radiation is generated by Bloch 
oscillations driven by the electric field £ applied in par- 
allel with the growth direction, i.e. with the superlat- 
tice lattice vector. Under influence of the electric field 
the Wannier-Stark ladder of quasi-stationary states is 
formed. The energy of emitted photons is determined 
by the separation between neighboring levels of the lad- 
der, hujBo = \e\£a, where lobo is the Bloch oscillation 
frequency and a denotes the period of the superlattice. 
Even though the idea of Bloch oscillations is old, it took 
a long time to find their experimental evidence i^i^ii^'^ 
For brief review see, e.g., Hartmann et ali^ and Leo4 

The strong transversal magnetic field has been used to 
quantize the free in-plane electron motion and to convert 
the quasi-three-dimensional electron structure of a super- 
lattice to quasi-one-dimensional Landau subbands, with 
the aim to change the electron dynamics and improve 
the condition for the terahertz emission, as described by 
Patane et al.,^ Scalari et al.,^^ and references therein. 

Completely different approach was used in our recent 
publication,— where we have theoretically studied the 
influence of the strong in-plane magnetic field on the 
electronic structure of superlattices. We have considered 
the superlattice subject to crossed electric and magnetic 
fields, £ and B, both applied "in-plane" , i.e. perpen- 
dicularly to the modulation direction, as an alternative 
source of radiation. In that case electrons are driven in 
parallel with the lattice vector by the Lorentz force. As 
electrons tunnel through the barriers, the cyclic motion 
along the electric field direction is superimposed to their 
otherwise straight-line drift due to the Lorentz force. The 
corresponding terahertz frequency uj^^^ is related to this 
cyclic motion and depends not only on the electric field 



£, as in the case of Bloch oscillations, but also on the 
applied magnetic field B. The frequency is given by 
WBii = 2TTVd/a where Vd = £ / B \s the electron drift ve- 
locity. 

Terahertz oscillations in still another configuration of 
the crossed fields was investigated experimentally by 
Qureshi,^^ who employed the transversal electric field, 
as in the standard Bloch configuration, combined with 
a strong in-plane magnetic field. The chaotic dynam- 
ics of electrons in the presence of the crossed fields and 
in the tilted magnetic fields was theoretically studied in 
papers li^ii^ii^ 

For Bloch oscillations, the very important issue is the 
intersubband tunneling of the Zener type. Lot of efforts 
were spent to clarify the condition under which Bloch 
oscillations could be observed and this problem is not 
yet definitely solvedJ^iiiiiSiii^ The same question needs 
to be addressed to the alternative magnetic-field-induced 
oscillations as we did not pay enough attention to this 
point in our previous publication>ii 

We will study the electrons in the presence of crossed 
magnetic and electric fields and subject to a unidirec- 




FIG. 1: The orientations of the electric and magnetic field 
with respect to the superlattice structure. Configurations (a) 
and (b) are considered in sections III. and IV., respectively. 
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tional potential with the period determined by the lat- 
tice vector a oriented in the growth direction of the su- 
perlattice. The considered geometrical arrangements of 
the fields and the superlattice are shown in Fig. [TJ The 
magnetic field acts in the direction z with a constant in- 
tensity B and the electric field of a constant intensity £ 
is parallel with the y axis. Two directions of the lattice 
vector are considered with respect to the orientation of £: 
the parallel one, a = (0, a, 0), with the periodic potential 
potential V{y) = V{y + a), and the perpendicular one, 
a = (a, 0, 0), with the periodic potential V{x) = V{x+a). 
The former configuration corresponds to the Bloch oscil- 
lations, the latter one to the alternative oscillations. 

While the three-dimensional superlattices subject to 
in-plane electric and magnetic fields have been proposed 
as an alternative source of the terahertz oscillations, 
in our theoretical analysis the simpler two-dimensional 
model is employed. As the z-dependent part of the three- 
dimensional Hamiltonian, pi /2m, does not play any role 
in the following theoretical consideration, we can omit it 
for simplicity without loss of generality. 



II. A FREE ELECTRON IN CROSSED FIELDS 

We start with description of a free electron in crossed 
fieldsi^ In general, the two-dimensional Hamiltonian of 
an electron in the electromagnetic field reads 



where 



1 / -\2 

H ^ \p- eA] +€( 

2m \ I 



£ = — grad0, i3 = curlj4, 



(1) 
(2) 

= —£ y and 



and e — — |e|. If we further assume 
employ the Landau gauge of the vector potential A 
(— Sy, 0,0), the Hamiltonian ([T]) reduces to 
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H^ = — (j>.~V\Byr + ^^-Wy 



2m 



2m 



(3) 



Since Hq commutes with p^., the separation of variables 
in the corresponding Schrodinger equation is possible. 

The eigenvalues and eigenstates of Ho are only slightly 
different from the more frequently described zero-electric- 
field case 



2m 2m 
The eigenfunctions of Hq can be written in the form 

1 



V'Onfc, (r) = 



(4) 



(5) 



where ^ j 2'k/ L^, j = 1,2, ■ ■ ■ , N and the eigenener- 
gies read 

Eon{k,)^hLu(n+l-] , w= n = 0,l,---. (6) 

\ 2 I m 



Lx and Ly denote the sample dimensions. 

The functions fonk^ {y) are the eigenfunctions of the 
harmonic oscillator, 



Vank^ (y) 



1 



^/2"nly/^e 



exp 



2P- 



Hr, 



y-yo 



(7) 

where Hn are the Hermit polynomials and £^ = h/mco. 
The coordinate yo = (-"^kx relates the center of the cy- 
clotron orbit, given by {^ponk^lyl^onkj = J/o, to the wave 
vector kx, yo € (0, Ly). Since k^ G (0, Ly/i"^) the number 
of states in a Landau level is given by A^ = LxLy/2Tr£'^ 
and the level degeneracy reads 



A^ 



\e\B 
h ■ 



(8) 



Very similar form of the eigenstates of the full Hamil- 
tonian Hq is obtained if we sum in Eq. ^ the terms 
linear and quadratic in y. Then Hq takes the form 



1 2 

Hq = — (px - mvd - \e\By) 
Zm 

pI , . mv'i 

+ ^ + Vd{Px - mvd) + -Y-, 



(9) 



in which the drift velocity Vd = £/B of an electron in 
crossed fields is employed. 

It follows from the similarity between expressions 
dH) and © that the eigenenergies of Hq can be con- 
structed from the zero-electric-field ones, EQn{kx), given 
by Eq. ([6]). They can be expressed as 



Eon{kx) = Eon{kx) + hkxVd + 



(10) 



where kx stands for kx — mvd/h. The Landau levels 
become tilted and, therefore, their degeneracy in kx is 
lifted. The eigenfunctions fQnk iv) remain practically 
unchanged, only the wave vector kx was replaced by kx 
in Eq. ©. 

In comparison with the zero-electric-field case the diag- 
onal matrix element of the velocity component is nonzero 
and given by 



iVonkJVxl'PQnkJ = T 



IdEjkxl 
h dkx 



Vd, 



(11) 



and the center of mass of the cyclotron orbit reads 

i^Onk^ \y\^0nkj = (^kx = 2/0 - Vd/cJ. 

This well known example demonstrates that the elec- 
tric field does change the zero-electric-field electronic 
structure only in an unessential way. Both systems are 
described by almost equivalent Hamiltonians which differ 
only by the shift of the coordinate center and the changes 
of the constants. We will show in the following sections 
that it is valid also in the presence of the superlattice 
periodic potential. 
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III. THE ELECTRIC FIELD PARALLEL WITH 
THE LATTICE VECTOR 

In the zero-magnetic-field limit this geometrical ar- 
rangement would lead to the standard Bloch oscillations. 
In crossed electric and magnetic fields the Hamiltonian 
can be written as a sum of Hq, as given by Eq. and 
of the periodic potential V{y), 



H = Ho + V{y). 
The resulting Hamiltonian reads 



(12) 



H 



2m 
2m 



mvd 



e\BvY 
V[y) + Vd{px - mvd) 



(13) 



Similarly as in the case of a free electron described above, 
it is obvious that the eigenstates of H are closely related 
to the eigenstates of the zero-electric-field Hamiltonian 



2m 



\e\Byf + ^ + V{y). (14) 



2 m 



For a given both Hamiltonians reduce to one- 
dimensional ones, which differ only by a constant Vdkx + 
mv'^/2 and a shift mvd/h, in the wave vector kx- There- 
fore both Hamiltonians yield essentially the same spec- 
trum of eigenvalues. 

The eigenenergies of the Hamiltonian (fT3| can be writ- 
ten, in analogy with Eq. pH)) . as 



Es{kx) = Esikx) + fikxVd 



(15) 



It follows from Eq. ([TS)) that Esikx) is a step-like func- 
tion with the step length K and the step height |e|£a, 
i.e. we can write 



Es{kx + K) = Es{kx) + hL0BO, t^BO 



|e|£o 

IT' 



(19) 



where lobq denotes the frequency of Bloch oscillations. 

This step function replaces the standard Wannier- 
Stark ladder obtained in the zero-magnetic-ficld case. 
Note that Esikx) are the true eigenenergies and not the 
resonances as the Wannier-Stark states, i.e. the magnetic 
field completely suppresses the Zener tunneling. 

The energy spectra of the above Hamiltonians can be 
found either by the direct numerical solution of the corre- 
sponding one-dimensional Schrodinger equations, or we 
can look for the eigenfunctions in the form of the linear 
combination of the free electron functions ([7]). As the 
resulting eigenenergies are the functions of fc^, we can 
evaluate the matrix elements for a given kx and diago- 
nalize the resulting matrix. 

The matrix elements of the Hamiltonian Eq. (I14p . 
are diagonal in kx- The periodic potential V{y) can be 
expanded into the Fourier series, 



V{y) = Y.^,cos{G,y), 



(20) 



with Gj = jG, G = 27r/a being the reciprocal lattice 
vector. Therefore, we can look for the eigenfunctions 
^sfcx (y) in the form 



~ ^ ^ Cn' IpOn' ; 



(21) 



Note that for the zero-electric-field case the periodic 
potential Vly) removes the degeneracy of Landau levels. 
The eigenenergies of the Hamiltonian become the 
Landau subbands Es{kx), s = 0, 1, 2, • • •, periodic in kx 
with the period 



K 



a 



a\e\B 



(16) 



We replace n by the integer s to distinguish between a 
subband and a level. The corresponding eigenfuctions 
ipsikx) are, similarly as iponk^ , bounded in y direction for 
any kx due to the confinement by the parabolic potential 
generated by B. 

The diagonal matrix elements of the velocity compo- 
nent Vx in crossed fields are now given by 

,^ I 1^ , IdEjkx) I dEsjkx) , ,^„. 
and the center of mass of the cyclotron orbit reads 

ivs-k. \y\^skj -yo- vd/uj - (18) 



and the matrix elements of H can be written as 
Hk 



huJ [n + -] 5n,n' + Tfcx 



(22) 



where Tk^^n,n' is the matrix element of a component of 
the periodic potential 



Tk^ 



1 



\/2"2"'n!n'!7r£2 





{y-ya?' 


j cxp 





y-Hn ( I f ^^-^ ) cos (G,y) dy. (23) 



Introducing the dimensionless variable rj = {y — ya)/(.^ we 
obtain 



V,- , = / exp I 

^ \/2"2"'n!7i'!7r 7-oo 



xH„ir,)Hn' iv) cos[G,h] + Gj2/o)d?7. (24) 
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To reduce this expression to a pair of standard tabular 
integrals we can employ the simple trigonometric relation 

cos{Gjh'i + GjUo) = 

cos{Gj£r]) cos(Gj2/o) ~ sm{Gj£ri) sin(Gjyo)- (25) 

After substitution, the examined T}.^^n,n' is divided into 
two parts Ii^n,n' and l2.n,n'- The first part, which in- 
cludes the cosine in the integrand, is equal to zero if the 
integrand is the odd function. For the even integrand it 
reads 



Il,n,n' 



^ _2cos(Gjyo) 
4^ ^ \/2"2"'n!n'!7r Jo 



exp 1—1] 



xHn{T])Hn'{v) cos{Gj£ri)dr]. 



(26) 



The explicit analytic form of Ii „ „/ (see,e.g.)2I) is given 
by 



h. 



2"n! 



■(G/)" -"exp 



G 



21)2 



The similarity between Hamiltonians p2l) and ([30)1 con- 
firms the validity of the expression p^ . 

As a simple example^^ we present the case of a weak 
perturbation V{y) = V cos(Gj/) which does not mix the 
Landau levels with different n. Then, replacing n by s to 
stress the difference between the level and the subband, 
the eigenenergies for the zero-electric-field Hamiltonian 
become 



Esikx) = huj [s+^ 



exp 



G2£2 



G^f^ 



cos I ^ ) , (31) 



and, finally, we get in accord with Eq. (fT5 



— / 1 \ — TflV 

Esikx) = Tiuj [s + - \ + hkxVd + 

( (G^e\ f2nk 

-h^exp I — I Ls ( ^— I cos 

for the eigenvalues in the electric field. 



j , (32) 



2/)2 



Gjf 



cos(Gjyo), 



(27) IV. THE ELECTRIC FIELD PERPENDICULAR 
TO THE LATTICE VECTOR 



for n'—n = 2p. Here ~" are the Laguerre polynomials, 
and p is an integer number. The second part l2,n,n' differs 
from the first one only by the sine written instead of the 
cosine 



l2,n,n' 



•^y^ 2sin(Gjyo) 



V2"2"'n!n'!7r Jo 

J 

xHn{r])Hn, (?/) sin(Gj^77)dr7. 



exp ( — 77 



(28) 



The part l2,n,n' is also given by the tabular integrals 
which are nonzero for the odd integrands 



l2,n,n' — 



I 2"n! 
2"'n'!' 



-1)- 



■{G,£) 



n'— n — 1 



2/72' 



exp 



Gji 



G2£2 



sin(Gjyo), 



(29) 



for n' — n ~ 2p + 1. Note that kx enters the matrix 
elements Tk^^n,n' through the identity Gjj/o — j'^T^k^/K. 
Only the minor changes are introduced by applying the 
electric field. The matrix elements of the Hamiltonian 
(fT3l) read 



1\ - mv 
Tiuj \ n + -\ + hkxjVd H — 



(30) 



This arrangement corresponds to the alternative oscil- 
lations. The Hamiltonian of an electron in crossed elec- 
tric and magnetic fields with the periodic modulation in 
the X direction can be written as a sum of the Hamilto- 
nian Hq, as given by Eq. and of the potential V{x), 



H = Ho + Vix), 



(33) 



or, explicitly. 



mvd 



\e\ByY 



+ 7; h V{x) + Vd{px - mvd) 

zm 



(34) 



We expect that due to the strong confinement by the 
magnetic field the eigenfunctions will be localized in y 
direction. For £ ^ 0, the eigenenergies should reduce 
to Es{ky), periodic in ky with the period K. This func- 
tion must be equivalent to Egikx) found in the previous 
section, as the result should not depend on the choice of 
calibration of the vector potential. 

Since the Hamiltonian p4p cannot be simply reduced 
to the one-dimensional one, we will look for the eigen- 
solutions starting from the linear combination of eigen- 
functions of Hq. The potential V{x) will be, similarly as 
in Sec. IHH considered in the form of the Fourier series 



V{x) 



3 



Vj cos{Gjx) 



(35) 



5 



The functions ips^^ (r), we are looking for, can be written 
as a Bloch sum, 

^sk^ (r) = C"'/e''=-'>o«'fc,,., (y), (36) 



where hj^j = + Gj, i.e. as the Hnear combination 
of functions the centers of which, I'^kxj, are arranged 
periodically along the y axis with the distance £"^0. 

For a given kx, the matrix elements of the Hamiltonian 
(l34l) can be written as 



(37) 



hk 



■xj' 



Here Tnj^n'j' denotes the matrix elements of the poten- 
tial V{x). They are products of the matrix elements of 
Vj cos{Gjx), calculated between the plane-wave parts of 
the wave functions, and of the overlap integrals of lo- 
calized functions ifon'k liv)- Oiily the matrix elements 
Vj-ji are nonzero. The centers of overlapping functions 
are ("^kxj and i'^kxj', respectively, and their distance is 



PGj-j,. Then T, 



nj,n' J' 



are equal to 



T I I —V- ■/ 



2"'n! 

r exp 

2"n'! ^ 



for n' > n and 



T ft —V- 



2"n'! 
2"'n! 



exp 



(38) 



LI] 



(39) 



for n' < n. They decrease exponentially with the dis- 
tance between the centers of functions fonk^. (y) and 

^On'k^^,iy)- 

Note that the electric field £ (the drift velocity Vd) and 
the wave vector kx enter only the diagonal part of the 
Hamiltonian, and, therefore, only the diagonal elements 
of the corresponding matrix equation which reads 



E 



0. 



(40) 



This allows us to predict an interesting kx dependence of 
the eigenstates of the above equation. 

First, it follows from the form of the matrix elements 
of the Hamiltonian p7|) that the eigenenergies depend 
linearly on kx^ and that the electrons move, in spite of the 
presence of the superlattice potential, with the same drift 
velocity Vd as free electrons in crossed fields. Moreover, if 



we replace kx hy kx + G, we can rewrite the Hamiltonian 
matrix elements p7p to the form 



TT_ 



kx -nj^nj 



Then equation 

E {^kx+G,nj,n'j' 
n'j' 

can be replaced by 



E). 



fiGvd- 



-n'j' 



n'j' 



(41) 
(42) 

(43) 



where E' stands for E—hGvd- Therefore we can conclude 
for the corresponding eigenvalues that 



Esikx + G) = Esikx) + hujB« , 



27rvd 



(44) 



The eigenvalues of Eq. P0)l should reduce for £ ^ to 
the zero- field eigenvalues Es periodic in ky. 

We will again illustrate the above consideration on the 
simple case of a weak perturbation y cos(G'a;). 

We start with the case of the^zero electric field^^ to 
show how the ky dependence of E^ is introduced. 

The function 

V-sfc. (r) = Cj'e'''="^-"=Vonfc,^., (y) = Y ^J' 1'^-^") (45) 
j' j' 

is assumed not to mix the Landau levels. Note that kx G 
(— 7r/a, ir/a) and j' 1, 2, • • • , jm where Jm = Lya/2'Kl'^. 
The equation 



{nj\H -E\^skx)^Q 



implies 



fiuj [ n + -] - E 



(46) 



(47) 



V 

-yexp 



I ) (Cj-1 +Cj + i) = 0. 



We write Cj in the form Cj = coexp(i(7j) and look for 
Co and q instead of Cj . The first consequence is that the 
solution of Eq. (|46p can be written in the form 



Es — huj 



f-j+T/cxp [-^) Ln[-^yosq. 

(48) 

It follows from (V's/c^ IV'sfc^) = 1 that JmCq = 1 and Cq = 
eG/Ly. 

If we write 'ipskx,q as 



we get from 



i'skx.q = co^e*«^'|nj'> 
i' 



{lpskx,q\^Pskx.,q') = 



(49) 



(50) 
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the expression q = {2t:v / Ly)£^G,y = 1,2,---. Taking 
into account that 2'kv / Ly = ky and Xq = i'^ky, we arrive 
at q — xqG — 2'Kky/K and finally 



Es{ky) = tiuj[s + l- 



+y exp 



cos 



(51) 



As expected, the change of the direction of the lattice 
vector has no influence and the only difference in com- 
parison with the V{y) case, see Eq. (|3T|) . is replacement 
of kx by ky. 

It is not too difficult to introduce the electric field into 
this simple model. First, the kx must be replaced by kx 
and q hy ky — Kq/2'iT in ipsk^,q- Then we can calcu- 
late the matrix elements of the Hamiltonian (|37p. which 
include the kx and Vd dependent parts. We get 



dk h' 



(52) 



and thus, taking into account the periodicity in fcj,, the 
eigenenergies of the corresponding equation can be easily 
calculated: 

1 f^- 

Er,S,k^ = Esiky)dky 

2 

+nkxVd+'^+rnujB^^, r = 0,±l,---. (53) 

In this example we limited our consideration to one sub- 
band Es{ky). Note that for more subbands the elec- 
tric field cannot cause the intersubband transition, as 
we use the electric-field-dependent basis of functions 
e*^-3'^(pp„,^^ ^ (y). Therefore, the solutions of Eq. (gD]) 
are the eigenstates and not the resonances. 



V. CONCLUSIONS 

The presented theoretical analysis has implicitly as- 
sumed a single electron model for a Bloch electron in 
crossed electric and magnetic fields, £ _L moving in 
a perfect superlattice crystal with the lattice vector a. 
Two geometrical arrangements were considered. 

In the case a || f , which would correspond to the Bloch 
oscillations for B = the magnetic field converts the 
standard Wannier-Stark ladder of resonances to the step- 
like eigenenergies, which are functions of the wave vector 
/c^, perpendicular to both £ and B. The step length 
is a\e\B/fi. The step height Tilobo corresponds to the 
frequency of the Bloch oscillations ojbo = |e|£a/?i. 

The alternative magnetic-field-induced oscillations 
were suggested as a possible source of terahertz radia- 
tior^ii for the case a _L f . The electron motion is com- 
posed from oscillations along a and the drift due to the 
Lorentz force with the velocity Vd = £ /B in the direction 
of k± . The resulting eigenstates resemble those of a free 
electron in crossed fields. The eigenenergies Em depends 
linearly on the fcj^, Em cx hvdk±, their separation on the 
energy scale is huj^^^ . 

The above conclusions are valid for the optimal con- 
ditions for coherent Bloch and alternative oscillations, 
neglecting the importance of additional scattering which 
can lead to damping of oscillations. 
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